library(methods)
library(stats4)
library(splines)
library(VGAM)

delta <- 0.9
bound <- 2097152
size <- 20971520
x <- rzipf(size,bound,delta)
#x <- sort(x)
sink("freq.dat")
for( i in 1:size)
	cat(x[i],"\n")
sink()
